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Abstract. For given convex set K in the plane (not necessarily bounded), we can con- 
struct the curve C for which the visibility (aperture) angle of this set has the same, 
prescribed value. We give the implicit formula for C, discuss some issues concerning prac- 
tical computations of C and bring several simple examples, when the equal visibility angle 
curves can be effectively computed explicitly. We conclude with some remarks about pos- 
sible directions for further research in this area. Extensive use of Open Source software 
(Sage, Pylab, IPython,. . . ) is a key feature of this article. 

Key words and phrases. Aperture angle, planar convex sets, equal aperture angle. 

t-H Krivka konstantnych uhlov viditel'nosti 

pre niektore konvexne mnoziny v rovine 

Abstrakt. Pre dami konvexmi mnozinu K v rovine (moze byt aj neohranicena) mozeme 
■^-j- zostrojit krivku C, pre ktoru je uhol viditelnosti tejto mnoziny konstantny, teda nadobuda 

predpisanu hodnotu. V clanku uvadzame implicitny vzorec na vypocet C, zamyslame 
sa nad praktickymi sposobmi vypoctu tejto krivky a uvadzame niekol'ko jednoduchych 
prikladov, kedy sa daju krivky rovnakych uhlov viditelnosti urcit explicitne. Nakoniec je 
naznacene, ake sii dalsie mozne smery vyskumu v tejto oblasti. V clanku je v podstatnej 
miere pouzity otvoreny softver (Sage, Pylab, IPython,. . . ). 

Kl'ucove slova. Uhol viditelnosti, rovinne konvexne mnoziny, konstantny uhol viditelnosti. 
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03 1. Introductory remarks 

Definition 1. Let Q be the convex set in the plane. The aperture (or visibility) 
angle <p(X, Q) of a point X Q with respect to the set Q is the angle of the 
smallest cone with apex X that contains Q (see e.g., [I], [2]). In this paper we 
will assume that the set Q is closed, but not necessarily bounded. 

For given angle a, < a < tt and the convex planar set Q we will denote 
by C(a, Q) the set of all points X for which the aperture angle <p(X, Q) has the 
value a. 

Example 1. Given the wedge W with internal angle for fixed angle a, three 
cases can occur: 

1. for a < the set C(a, W) is empty, 
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2. if a = <& then C(a,W) is the wedge V (W rotated by angle ir about its 
apex) , see Figure [T] 

3. for a > $ the set C(a, W) is the angle V with value 2a — opposite to 
wedge and having common bisector with angle i?, see Figure [2} 




Figure 1. C(a, W) = V for a = 0. Figure 2. C(a, W) = V for a > 

In the Example [T] the convex set W was unbounded and not strictly convex. 
The set C(a,W) was two-dimensional in one "singular" case. Next, we will 
consider unbounded, strictly convex set. 

Example 2. Let us define the convex set Q as Q = {(x,y) 6l 2 : y > x 2 }, 
i.e., Q is the parabola y = x 2 together with its interior points. 

In this case, it is not difficult to compute the set C(a, Q) explicitly. For 
arbitrary exterior point A = (x±,yi) Q we can construct exactly two tangent 
lines to Q passing through A. The slopes k\, hi of that tangent lines are 

ki = tan Lpi = 2 (xx + \Jx\- y^j k 2 = tan (p 2 = 2 \ x\ - \Jx\- y^j , (1) 

where tpi,<p2 are the slope angles of two tangent lines. For given angle a we seek 
the set of all points A such that 

tan - tan (/? 2 h - k 2 

tan(<£i - <f 2 ) = —— : = : , ; r = tan a = K. 2 

1 + tan (pi tan cp 2 1 + k\k 2 

If k\k 2 —1, we can substitute for fei, k 2 from ([T| and the solutions of resulting 
quadratic equation give us the explicit formula for the curve C(a, Q) 

. . -K 2 - 2 ± 2V4K 2 x 2 + ^2 + 1 

yp( x ) = » ( 3 ) 

where the plus sign is valid for angles a > tt/2 and the sign minus for sharp 
angles a. For a = tt/4 this formula is simply y p (x) = —(2\/4x 2 + 2 + 3)/4 and 
this is the curve with asymptotes y = ±x. For arbitrary K ^ 0, the asymptotes 
of curve given by formula Q are y = dzx/K. 

For the singular case k\k 2 = —1, i.e., for a = tt/2 we obtain immediately 
y p (x) = —1/4, which is the directrix of parabola y = x 2 . From the formula Q3J) 



M. Kaukic: Equal aperture angles for some convex plane sets 



3 



we can see that for angle a < tt/2 the function y p (x) is concave and for a > ir/2 
we obtain convex function y p {x). 

On the Figure [3] we have plotted the graphs of the functions y p (x) for several 
angles a. From the bottom to top curve, the angles are: 7r/4, 7r/3.5, 7r/3, it/2.5, 
it/2, vr/1.5, it/1.2. 

Now, the set C{a,Q) is nonempty 
for all angles < a < it. For compact 
(in M 2 - closed and bounded) Q the 
set C(a,Q) (if nonempty) forms the 
closed curve, which we will call con- 
stant visibility angle curve. This 
curve have been investigated mainly 
for the special case of convex polyg- 
onal sets Q (see [6] , [l] , |] ) . The 
papers cited are concerned with op- 
timization problems (maximization) 
for visibility angle. In the paper [3] 
the authors solved three-dimensional 

problem of maximization of visibility Figure 3 - Functions Wp(*) for parabolic region, 
angle for convex polyhedron, viewed 

from given line segment. Good sources of information about related problems are 
the Pirzadeh's Master thesis Q and the article jjjj. 




2. Constant visibility angle curves for convex functions 

We can, in principle, compute the visibility angle curve for arbitrary convex 
function. Let us assume that the function y = f{x) is strictly convex on finite 
segment I = [a, b). Further, we assume that f(x) has continuous first and second 
derivative on /. 

Given the (arbitrary, but fixed) point x G I and the angle a, the problem is to 
find another point x 6 I such that (cf. equation Q) 

^ ~ ™ — tan a = if or /'(,) = = M. (4) 



1 + f'(x)f(x) 



1 + Kf(x) 



For given x, the constant M is known, therefore we have to solve simple equation 
F{x) = f'{x) — M = 0. For some simple functions we can do that analytically 
but, in general, it is necessary to use one of suitable numerical methods. For 
example, using Newton method we get the iterational sequence: 



choose xq € I, x ra +i 



/'(an) - M 

f"(Xn) 



for n = 0, 1, . . . 
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We will not analyze the convergence assumptions here. It will be the subject of 
further research. In the next section we give two nontrivial examples of computing 
the constant visibility curve for smooth, bounded convex sets. 

3. Constant visibility angle curves for sine and ellipse 

Example 3. Let Q be the convex region, described by 

Q = {{x, y) e R 2 : < x < tt, - sin x < y < 0}. 

We will take into account only visibility "from below", i.e., from points with 
negative y coordinate. 

On the Figure [4] we can see four main 
regions in the lower halplane. In each of 
them the computation of C(a, Q) is a bit 
different. From the symmetry of Q it fol- 
lows that it is sufficient to investigate only 
parts of regions with x < tt/2. In the re- 
gion (I) there are only the points of curves 
C(a, Q) for a > tt/2. In the regions (III) 
resp. (IV) the visibility angles curve al- 
ways contain points A resp. B. For the 
points in region (II) we have a < ir/2, the 
set Q is visible only as the segment AB, 
Figure 4. Regions for computing C(a, Q). tmis t h e curves C(a, Q) are circular arcs. 

In this example, the curves C(a, Q) can be computed analytically, but to com- 
pute them "by hand" as in Example ^ is very tedious and error-prone. For 
symbolic computations, we utilized Open Source system Sage [8]. For numeri- 
cal experiments the programming language Python j9], with user friendly shell 
IPython [12] and modules Numpy [TO], Scipy [7] were used. Finally, nearly all 
graphics in this article was generated by Python module Matplotlib |11| . 

Programs used for computations and generation of graphics can be downloaded 
from [http : // f eelmath. inf o/images/constangle/constang le .zip| For com- 
putation of points of the C(a, Q) curves (except circular arcs) we need: 

1. take an arbitrary point £, < ^ < tt/2 and the tangent line t\ to Q in the 
point Ti = (£,-sin£), 

2. compute the associated tangent line £2 to Q such that Z(ii,t2) = a, 

3. the intersection point T = (x(£,),y(l;)) of tangent lines t±,t2 lies on the 
curve C(a, Q). 

The behaviour of curve C(a,Q) depends on the angle a and it is very inter- 
esting. For a near (but less than) straight angle, those curves are convex. For 
certain obtuse angles close to tt/2 the central part of curves (i.e., the part con- 
tained in the region (I)) have waveform shape. Minimal amplitudes of "wave 
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oscillations" we observed for the angle of about 1.9 radians. For obtuse an- 
gles, smaller than (approximately) 1.88 radians the above-mentioned central 
part is concave. On the figure [5] we can see the curves C(a,Q) for ten an- 




Figure 5. C(a,Q) for obtuse angles. Figure 6. Wawe-shaped central parts. 

gles aj~ = 7r/(l + 0.1 k), k = 1, 2, ... , 10. Figure [6] shows central waveform parts 
of twenty curves C(a,Q) for angles in interval from 1.92 to 1.95 radians (the 
picture is zoomed vertically). For sharp angles the curves C(a,Q) are not very 
interesting, for a — > they are close to big circular arcs with central angle 2a. 

Example 4. We take the convex region Q, bounded by ellipse 

x = a cos ip, y = b sin ip, < (p < 2ir, a > b. 

Here, using the Sage system, we can derive the formula for C(a,Q) in polar 
coordinates. Note, that the curve C(a,Q) for a = tt/2 is simply the circle with 




Figure 7. C(a, Q) for "thin" ellipses. Figure 8. C(a, Q) for "thick" ellipses. 

radius r = Va? + b 2 . Obviously, for angles close to straight angle, the constant 
view angle curves are convex. If the ratio b/a is greater than approximately 0.707 
(hypothesis: maybe this value is y/2/2), then the regions bounded by curves 
C(a,Q) are convex for all a. On the other side, if the ratio b/a is small, the 
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curves C{a,Q) for small angles will approach the two circular arcs with central 
angle 2a. 

Directions for further research 

In this paper, we brought only simple examples, but already we have seen some 
interesting properties of curves of equal aperture angle. The broad research area 
is to find efficient numerical algorithms for computations of those curves. It 
would be interesting to approximate C(a,Q) with the corresponding curve for 
polygonal approximation of boundary of region Q. Here we can explore some 
results from [l] or [6] (e.g., the method of rotating wedges). 
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